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Abstract 

I I We study the integer minimization of a quasiconvex polynomial with 

quasiconvex polynomial constraints. We propose a new algorithm that is 
an improvement upon the best known algorithm due to Heinz {Journal 
• of Complexity, 2005). This improvement is achieved by applying a new 

modern Lenstra-type algorithm, finding optimal ellipsoid roundings, and 
considering sparse encodings of polynomials. For the bounded case, our 
g algorithm attains a time-complexity of s(r/Md)°(i)22"i°S2(n)+o(») ^hen 

I I M is a bound on the number of monomials in each polynomial and r is the 

binary encoding length of a bound on the feasible region. In the general 

m case, s/0(i)^o{n)22«iog2{n)+o(n)_ assume d > 2 is a bound on 

^ the total degree of the polynomials and I bounds the maximum binary 

encoding size of the input. 

1 Introduction 

We study the integer minimization of a quasiconvex polynomial with quasicon- 
vex polynomial constraints. That is, given F,Fi, . . . ,Fg G Z[x] = Z[xi, . . . , a;„] 

. . quasiconvex polynomials with integer coefficients, we wish to solve the following 

^ problem 

k> min -f'(x) 

}_J subject to Ei(x) < for alH = 1,. .. ,s (1) 

<^ X e Z". 

A function F : M" ^ M is called quasiconvex if for every a S M, the lower level 
set {x e M" : F{x) < a} is a convex subset of R". Some quasiconvex programs 
reduce nicely to convex programs, see for instance, [9], but this is not likely to 
be the case in general. Studying quasiconvex integer minimization opens up a 
larger class of functions that we can optimize over. 

We approach the optimization problem by setting Fq = F — z* and solving 
the feasibility problem over F n Z" , where 



y := {x e M" : Fiix) < 0, i = 0, 1, . . . , s} , (2) 
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and applying binary search on objective values until we find an optimal solution. 
Strict inequalities are used to ensure that if Y is non-empty, then it is full 
dimensional in M". Since -Fi(x) e Z for all x e Z", problem (1) can be easily 
formulated by weak inequalities. This follows from the observation that the 
inequalities z <0 and z + 1 <0 are equivalent for 2; G Z. 

We use a modern Lenstra-type algorithm to solve the integer feasibility 
problem. Lcnstra's algorithm was the first algorithm to solve integer linear 
optimization in polynomial time when the dimension is fixed. It can be applied 
to any family of convex sets C in M" provided that we can solve the ellipsoid 
rounding problem over sets in C. Khachiyan and Porkolab [23] showed that 
Lenstra's algorithm could be generalized to operate on convex semialgebraic 
sets, having time-complexity of /'^^^^(srf)'^^" \ For the specific case of quasicon- 
vcx polynomial minimization, the current best algorithm is due to Heinz and 
has time-complexity of ,si*^^^-'d*^'-"^2'^(" \ where d > 2 is an upper bound on 
the total degree of the polynomials and I is the maximum binary encoding size 
of all coefficients. 

Our improvement over Heinz's algorithm comes primarily from the modern 
Lenstra-type algorithm that we present. Heinz developed a shallow cut separa- 
tion oracle to show that Lenstra's original algorithm applies to the quasiconvcx 
minimization problem (1). We generalize Heinz's shallow cut separation oracle 
to show that the modern Lenstra-type algorithm works for the quasiconvex min- 
imization problem (1). Wc also provide a structure of evaluating polynomials 
that exploits sparsity, which allows us to state a more precise complexity of the 
algorithm based on the number of monomials given in the input. 

Theorem 1.1. Let F, Fi, . . . ,Fs G Z[x] be sparsely encoded quasiconvex polyno- 
mials. Let d > 2 be an upper bound for the degree of the polynomials Fq, . . . , Fs, 
let M be the maximum number of monomials in each, and let the binary length 
of the coefficients be bounded by I. Then there exists an algorithm for the min- 
imization problem (1) which computes a minimum point or confirms that such 
a point does not exist. 

(a) // the continuous relaxation of the feasible region is bounded such that r 
is the binary encoding length of a bound on that region with r < ld^^"\ 
then the algorithm has time- complexity 0/ s(r/Mrf)'^(^^22"'°S2(")+0(") and 
output- complexity of {I + r){dn)^^^\ 

(b) Otherwise, the algorithm has time- complexity 0/ sZ'^(-^)cZ'^(")2^"^°S2(") and 
output- complexity of ld^^'^\ 

For d = 0(1), this complexity is 5/0(1)22" i°g2(»)+o(n)_ 

If = 0(n'=) for some A; > 0, then the complexity becomes sZO(i)2°("'°s(")). 

Lenstra's algorithm solves the integer feasibility problem for a convex set Y 
by first finding a pair of concentric ellipsoids, E,E' = j^E such that E' CY C E, 
where E' is a scaled version of E with respect to the center. If wc can determine 
that E' n Z" is non-empty, then we are done. Otherwise, we find a direction 
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Figure 1: Main idea of Lenstra's Algorithm 



of minimal width and branch into integer hyperplanes that cover i?, creating 
lower dimensional subproblems to solve. The same approach is then applied to 
solving each lower dimensional subproblem. 

The complexity of Lenstra's algorithm is guided primarily by the number of 
subcases that it must evaluate. As shown in section 4, the number of subcases 
in each step is bounded by 2u! {n)f3{n) + 3 where uj{n) is called the lattice width 
direction of the inner ellipsoid and /3(n) is the rounding radius or scaling between 
the concentric ellipsoids that are obtained. The lattice width uj{n) will be used 
to determine if E' n Z" ^ 0. This leads to a total number of subcases bounded 

by 

n 

total # of subcases < 2uj{i)j3{i) + 3. 

1=1 

The idea of the modern Lenstra-type algorithm is not new. Kannan was the 
first to reduce the number of subcases by solving lattice problems. Recently, 
extraordinary new lattice algorithms and bounds from the geometry of numbers 
reveal a better complexity. A new idea we present is finding an ellipsoid rounding 
with an optimal rounding radius of Pin) = 0{n). This is done by using results 
for the maximal radius of a ball inscribed in the convex hull of m points on the 
sphere of radius 1. We will explain these improvements and how they affect the 
complexity of our algorithm. 

In section 2 we will discuss ellipsoid rounding from the point of view of the 
shallow cut ellipsoid method and show how we improve from Oirt'/'^) to 
0(n), which has never been done before in Lenstra's algorithm. This improve- 
ment allows for a better exponential coefficient in the resulting complexity. 
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In section 3 wc will explain how Kannan's improvement of Lenstra's al- 
gorithm reduces the number of subcases exponentially. This section will begin 
with a discussion of lattice theory, flatness directions, and the geometry of num- 
bers, and wc will reveal that u}{n) can be improved from 2*^*^" ^ as used in the 
original Lenstra algorithm [28] along with Khachiyan and Porkolab [23] and in 
Heinz [15], to 0{n) as stated in Eisenbrand [12]. The important feature that 
wc point out is how new lattice algorithms allow this computation to be done 
determistically in 2'^("^ time, as opposed to 2'^("''°s"), creating a better overall 
complexity. 

In section 4 wc state our modern Lenstra-type algorithm more precisely and 
explain its complexity. 

In section 5 we will discuss polynomial encoding and quasiconvex polyno- 
mials and then show how to make shallow cuts for sets given by quasiconvex 
polynomials. This allows Lenstra's algorithm to be applied to our problem. 

In section 6 we discuss the proof of Theorem 1.1. 

2 Ellipsoid Rounding 

The flrst step of Lenstra's algorithm is to flnd a pair of concentric ellipsoids, 
one inside and one containing the feasible region. We will write ellipsoids in the 
form E{A,a.) = {x € M" : ||x-a||A-i < 1} where ||v||b := VWBv, A e M"><" 
is a positive definite matrix and a e M". For example, E{a^I,0) = B{a,0) is 
the ball of radius a centered at the origin. Let K be a convex set. The ellipsoid 
E{A, a) is a (3-rounding of Y if 

E{j,A,ai) CY CE{A,Si) 

where /3 is called the radius of the rounding [32]. John [18] showed there exists 
a n-rounding for any convex set. Conversely, for a simplex, a n-rounding is the 
best possible rounding. Finding an optimal rounding will reduce the number of 
subcases that need to be analyzed in Lenstra's algorithm. In section 4 we will 
show precisely where /3 affects the complexity of integer optimization. In this 
section we will explain how to construct a 0(n)-rounding. 

There arc several methods for ellipsoid rounding. Nesterov describes an 
algorithm to obtain a 7n-rounding, 7 > 1 for an arbitrary convex set and also 
how to obtain a 7-;/n-rounding for centrally symmetric convex sets [32], but 
each is based on the assumption that a difficult optimization problem can be 
solved. For a specific case, Nesterov uses linear programming, whereas we would 
need to maximize over nonlinear polynomials. In our model, no supplementary 
optimization problem need be solved. Ellipsoid roundings have also been studied 
by Khachiyan [22], which was improved by [27] and [35]. Some other methods 
use a volumetric barrier [2,3]. Unfortunately, none of these have been shown to 
round general convex sets. 

We will use the original approach, which is to employ the shallow cut ellipsoid 
method [13]. This can be applied to any class of convex sets C for which there 
exists a shallow cut separation oracle. 
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Definition 2.1 (Shallow cut separation oracle [13]). A shallow cut separation 
oracle for a convex set F C M" is an oracle which, for an input a S Q" and a 
rational positive definite matrix A, outputs one of the following: 

1. verification that E{A,a.) is a /3-rounding of Y or, 

2. a vector c e Q" such that the half-space 

|x e K" : c'^x < c'^a + ^^^Vc^A-icj D F n E{A, a). 

Tiieorem 2.2 (Shallow Cut Ellipsoid Method [13]). There exists an oracle- 
polynomial time algorithm that for any number e > and for any circumscribed 
closed convex set Y C B{R, 0) given a shallow cut separation oracle finds a 
positive definite matrix A e Q"^" and a point a G Q" such that one of the 
following holds: 

(i) E{A,a) is a ji-rounding ofY. 

(ii) Y C E{A,Si) and voI{E{A,sl)) < e 

This algorithm runs in time oracle-polynomial in n+ {R) + (e) . 



The main difficulty in generalizing Lenstra's algorithm to different classes 
of convex sets is creating shallow cuts. Khachiyan and Porkolab show that for 
an arbitrary semialgebraic set, a shallow cut can be computed in l^^^\sd)'-'^'^ ^ 
time [23, Lemma 4.1]. In our algorithm, wc intend to do much better for the 
specific case of quasiconvex polynomials, although this complexity discussion 
will be saved for section 5 when we discuss quasiconvex polynomials. 

Suppose that we have an ellipsoid E{A, a) D Y. Ellipsoids are afSne trans- 
formations of the unit ball; that is, there exists an afhne map t : M" — > M" such 
that t(£'(A, a)) = B{1,0) D t{Y). The standard method of creating a shallow 
cut separation oracle is to observe that points in 5(^^,0) \ t{Y) will often 
admit a shallow cut. For the case when y is a polyhedron, finding such a point 
directly admits a shallow cut, whereas in the case of quasiconvex minimization, 
Heinz showed that with a little more work, we could find a shallow cut. This 
will be explained in detail in section 5. The remainder of this section will focus 
on realizing a /9-rounding. 

On the other hand, if we find a set of points V within the ball of radius 
that arc in t{Y), then any ball B{0, ^) C conv(y) will admit a rounding 

since r~^(B(-g,0)) = E{-^A,a.) CY c E{A,a). The rounding radius /3 is then 
dependent on the maximum inscribed sphere inside conv(y), as seen in Figure 2. 

Using the cross-polytope 

Vl = {±—^ei :i = l,...,n}, 
n+1 

where ej is the i*^ unit vector, [13] obtains a 0(n^/^ )-rounding of a polytope, 
which is also used in [23]. Heinz used this idea to obtain a rounding of a convex 
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Figure 2: Under the affine transformation, the open dots are test points for 
Heinz's separation oracle. This method avoids the numerical issues of finding 
exact square roots. 



region given by quasiconvex polynomials [15]. In order to overcome numerical 
issues of requiring exact arithmetic, Heinz chose the points 

Vh = {iAiG^ : i = 1, . . . , n} 

where ^_^3^^ < Ai < Heinz's choice also obtains an 0(n^/^)-rounding. We 
will generalize and improve Heinz's method by applying sphere approximating 
polytopes of Kochol [25] that attain an optimal bound within a constant factor. 
A note from Kochol, modified to give more detail, shows the following result. 

Theorem 2.3 (Theorem 3 in [26]). Let n, m be positive integers, 2n < m < c" , 
where c > 1 is a constant. Let p(n,m) as the maximal radius of a ball (with 
center at the origin) contained in the convex hull of m points chosen from the 
n-dimensional sphere of radius 1. Then there exist constants Ci and Ci such that 

/log(m/n) /log(m/n) 

Cl V < p[n,m) < C2\ . 

V n V n 

Furthermore, there exists a polynomial time algorithm in n and m to construct 
a set of vectors C Z" with \ V\ < m such that the polytope with extreme points 
{v/||v||2 : V e V} is symmetric across all axes and attains such bounds. 

Kochol notes that choosing m :— points improves the 0(ri'^/^)-rounding 
to 0{n^^^ /y/log n) and still allows a polynomial time rounding, and improves 
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Figure 4: An inscribed polytope and a rational approximation with vertices, 
respectively, as filled in dots and open dots. The innermost circle is the ball 
that we can guarantee will remain inside the convex hull of the approximated 
points. 
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upon the complexity for Lenstra's algorithm given in [13]. Theorem 2.3 demon- 
strates that using an exponential number of points is necessary to obtain an 
0(n)-rounding via the shallow cut ellipsoid method. A tighter rounding is ad- 
vantageous, even if it rcqiiircs an exponential mimbcr of evaluations. In our 
new algorithm we choose a single exponential number of points, m ;= n2" will 
sufBce, to obtain an 0(n)-rounding. This improves the exponential coefficient 
in the final complexity. 

We will now show that numerical approximations of Kochol's approximating 
spheres will still allow for an optimal rounding. We denote the sphere of radius 
r as ^""^(r) = {x e K": ||x||2 = r}. A l-net of ^""^(r) is a set of points 
N C S'"~^(r) such that for any point z e /S'"~^(r), there exists a point v such 
that ||v-z||2 < 1. 

Lemma 2.4. Let N he a l-net o/S'"-i(l), letO<a<\ and letQ<e<\-a. 
Suppose that N is an e-approximation of N, that is to say that for all v £ N 
there exists avGN such that 

||v - v||2 < e, 

then conv(iV) D ^"-^(a). 

Proof. Suppose there exists a point z € S"~^{a) that docs not belong to 
conv(iV), then separating z from conv(iV) by a hyperplane p^ we get a cap 
of which is disjoint from N and its top t where t is perpendicular 

to Since t is in S'"~^(l), there exists a point v E N from the l-net that 
satisfies ||t — v||2 < 1. See Figure 2 for the geometry. Letting d be the minimum 
distance between v and the hyperplane pz, we can see that d> g — a > e, which 
is a contradiction since N is an e-approximation of N. □ 

Corollary 2.5. Let N be the set of points given by Kochol's construction for 

an approximation of the unit sphere and let N be an e-approximation of N 
with < e < I . Then there exist a constant C\ (e) such that the ball of radius 
^^^J]og{^Jn) contained in conv(7V). 

Proof. Using Lemma 2.4 with a = j, the proof is very similar to Theorem 1 

in [25]. □ 

This is now the template for our separation oracle. We will choose m := n2" 
test points according to an approximation of Kochol's set of points. If all the 
test points are feasible, we obtain an 0(n)-rounding. Otherwise, we find an 
infeasible test point and generate a shallow cut. The specific algorithm for 
finding a shallow cut for quasiconvex polynomials will be presented in section 
5. 
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Figure 5: Geometry of the proof for Lemma 2.4 

3 Integer Feasibility and Subcases 

We assume now to have an elhpsoid rounding 

£;(^Aa)cy ci;(A,a). 

The next step in Lenstra's algorithm is to determine if the inner elhpsoid con- 
tains an integer point. A simple, yet powerful, way to do that is Khinchin's 
Flatness Theorem, which roughly states that if the minimum width of a convex 
body is greater than some constant a;(n), then the convex body contains an 
integer point. If the minimum width is less than uj{n), then we can branch into 
integer hyperplanes perpendicular to the minimum width direction. Since we 
must branch on the larger ellipsoid, we will have fewer than /3(n)aj(n) subcases 
to branch into. We will first review lattice theory and flatness directions, and 
present theorems for reducing the complexity of Lenstra's algorithm. 

3.1 Lattices 

Given m linearly independent vectors bi, b2, . . . , b^ e M", the lattice A gener- 
ated by i? = [bi, . . . , b„j] is the linear transformation of Z" 

A = C{B) {Bz : z e Z"}. 

The set of vectors bi, . . . , hm, or similarly B, is called a basis for the lattice A. 
Lattices naturally arise when looking for integer points in ellipsoids, since an 
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Figure 6: If the lattice width of the inner elhpsoid is greater than then we 

know it contains a lattice point. Otherwise, we are in the case displayed here, 
with roughly fewer than f3{n)uj{n) subcases to project into. We later use the 
loose bound of 2f3{n)uj{n) + 3 for the number of subcases, as shown later in this 
section. 

ellipsoid is an afBne transformation of -6(1, 0). We will need the following related 
properties of a lattice. The dual lattice A* is given by 

A* {v e span(B) : v^B e Z" }. 

In particular, one can show that when B is fuU rank. A* = C{{B-'^)'^). We will 
use a transference bound later, which is an inequality relating the properties of 
a lattice and its dual. 

The covering radius /i(A) is the smallest number a > such that the set of 
closed balls of radius a centered at each lattice point completely covers all of 
K". 

The shortest vector problem (SVP) is a well-studied and important problem 
in lattice theory with many applications. The goal is to find a non-zero lattice 
vector X e A \ {0} with minimal length A(A). For our purposes, we will find the 
shortest vector with respect to the Euclidean norm; therefore, 

A(A) = min llvlU = minlllvlU : v = Sz, z e Z" \ {0}} . 
veA\{o} 

Similarly, the closest vector problem (CVP) is to find the closest vector in 
the lattice to a given point t. Again, we will use the Euclidean norm for this 
problem. We will define CVP(A,t) as 

CVP(A,t) = argminllv-tlb. 

VGA 

For further review on lattices, see, for instance, [8] or [12]. 
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3.2 Lattice Widths and the Shortest Vector Problem 



Kannan first observed that SVP could be used to minimize the mimbcr of 
branching directions in Lenstra's algorithm [19]. We follow Eisenbrand in pre- 
senting this in the context of flatness directions [12]. Let K C M" be a non- 
empty closed subset of K" and let d e K" be a vector. The width of Y along d 
is the number 

WdiK) = maxjd^x : x e K} — min{d^x : x e K}. 
The lattice width of Y is defined as 

w(K) = min WdiK), 

dGZ"\{0} 

and any d that minimizes Wd{K) is called a, flatness direction of K. 

Theorem 3.1 (Khinchin's flatness theorem [24]). There exists a function uj{n), 
depending only on the dimension, such that if K cM." is convex and w{K) > 
uj{n), then K contains an integer point. 

The best known bound for w{n) is 0(n^/^) and it is conjectured that uj{n) = 
<d{n) [5]. We will see in the next subsection that, for the speciflc case of ellip- 
soids, we can obtain this bound. 

Flatness directions are invariant under dilations. This is easily shown for 
the case of ellipsoids. 

Lemma 3.2. Let d S Z" 6e a flatness direction for E{A,aL). Then for any 
13 gM, d is a flatness direction for E{j^A, a) with 

w{E{j,A,a)) = ^w{E{A,a)). 

Remark 3.3. For an ellipsoid, a flatness direction can be computed by solving 
the shortest vector problem in the lattice A = £((A^/^)^). To see this, consider 
the width along a direction d of the ellipsoid E{A, 0), 

Wd{E{A, a)) = max{d^x : x e E{A, 0)} - min{d'^x : x e E{A, 0)} 
= max d^fxi — X2). 

xi,X2e-E(A,0) 

We have d^(xi - X2) = d^A^/^{A-^/'^yii - A-^/^X2) where A-^/'^xi and 
A~^/^X2 are contained in the unit ball if and only if xi,X2 € E{A,0). Thus 
properly choosing xi and X2 on the boundary of E{A, 0), we see that 

Wd{E{A,0)) = 2||d^^V2||2 = 2\\{A'/Yd\\2. 

Finding the minimum lattice width is reduced to solving a SVP over the lattice 
A. After this transformation, integer points from £'(^A, a) fl Z" now live in 
S(l,^-V2a) n A* where A* = £(A-V2). 
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Figure 7: We use SVP on A = to solve for a flatness direction. 

If this shows that the lattice width w(£'(^j4,a)) > uj{n), then it contains a 
lattice point. We then recover a feasible integer point by solving CVP(A* = 
C{A^^/^), A^^^^a), where the target point is the center transformed by the 
lattice basis. The solution transformed back to the original space will be the 
closest lattice point to a with respect to the ellipsoidal norm, and therefore is 
contained in the ellipsoid. 



There are two computations arising here: 

1. We find a shortest vector in A to determine a flatness direction. If 2A(A) < 
u;(n), then we will project into a minimal number of subcases. 

2. If 2A(A) > uj{n) then we have confirmed that E{^A, a) contains an integer 

point. To recover this lattice point, we solve z — CVP(A*, A^^/^a) and then 
set X* — A^l'^-L. Since x* is then a closest integer point to a with respect to 
the ellipsoidal norm and we know that E{-^A, a) contains an integer point, we 
have X* e £:(^^,a)nZ". 

3.3 Results from the Geometry of Numbers 

The geometry of numbers produces a small bound on the lattice width of an el- 
lipsoid not containing an integer point. Using properties of LLL reduced bases, 
Lenstra originally observed that this value did not exceed 2*^^" ^ [28]. By con- 
sidering the specific case of ellipsoids, we can achieve an 0{n) bound. For an 
arbitrary lattice, the product of the length of a shortest vector in a lattice and 
the covering radius of the dual lattice is bounded by a constant f{n) dependent 
only on dimension. Using the Fourier transform applied to a probability mea- 
sure on a lattice, Banaszczyk showed that this function is bounded by a linear 
factor in the dimension n. 
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Theorem 3.4 (Theorem 2.2 in [4]). Let A C M" 6e a lattice with n > 1. Then 

A(A)^(A*) < /(n) < in. 

If we assume that a specific eUipsoid does not contain a lattice point, then 
the covering radius of the associated lattice is greater than one. Since the lattice 
width of an ellipsoid is simply twice the length of a shortest vector, we obtain 
the following inequality for ellipsoids. 

Theorem 3.5 (Theorem 14.26 in [12]). If E{A,a.) C K" ts an ellipsoid that 
does not contain an integer point, then w{E{A,a)) < 2f{n). 

A convenient bound follows directly from Theorems 3.4 and 3.5. 

Corollary 3.6. Let E{A,a.) c be an ellipsoid not containing an integer 

point, then w{E{A,a.)) < n. 

Hence, if E{-^A, a) does not contain an integer point, then w{E{A, a)) < n^. 

3.4 Complexity of the Shortest Vector Problem and the 
Closest Vector Problem 

The shortest vector problem (SVP) has been shown to be NP-hard, even to 
approximate within a constant factor [29]. Until recently, the best known 
deterministic solution to SVP was by Kannan, with time-complexity n! = 
20(niogn) |-2Qj r^Y^^ well-known Ajtai, Kumar, and Sivakumar [1] sieving method 
is a probabilistic method that solves SVP with very high probability and achieved 
the first single exponential time-complexity, which was shown by [33] to be 2^-^". 
Micciancio and Voulgaris improved this type of method to achieve a run time 
of 23 199" [31]. 

The closest vector problem (CVP) has also been shown to be NP-hard, 
even to approximate it within a polynomial factor [11]. Kannan presented an 
algorithm to solve CVP in n! time [20], and although there have been some 
improvements [14,16], none have been able to achieve a single exponential time- 
complexity. 

Micciancio and Voulgaris discovered the first deterministic single exponential 
algorithms for CVP and SVP [30]. This is an exciting and impressive result. 

Theorem 3.7 (Corollarys 4.3 and 4.4 in [30]). There are deterministic algo- 
rithms to solve SVP and CVP, with the Euclidean norm, that both have time- 
complexity n'^(-'^^2^". 

This result, for the first time, allows the complexity of SVP to be much 
smaller than the complexity of Lenstra's algorithm. 

3.5 Projection 

We will need the following simple lemma, which can, for instance, be found 
in [17]. We indicate it here with a proof to give a precise complexity. The proof 
uses the common fact that if i? is a basis for a lattice A, and U is a, unimodular 
matrix, then BU is also a basis for A. 
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Lemma 3.8. Suppose A C Z" is a lattice with basis {bi, . . . ,b„j} and m > 2. 
Suppose d S A \ {0} is primitive (i.e., ad ^ A for all < a < 1), and let 

Aj e Z for i = 1, . . . , m such that d = Aibi H h Ambm- Then there exists an 

algorithm that computes vectors h2, ■ ■ , b„j such that {d, b2, . . . , b„} is a basis 
for A. This algorithm has time- complexity (nlogL)'^^^^ where L is an upper 
bound on the absolute values of Xi and the entries of hi for all i = 1, . . . ,m. 

Proof. Let B = [bi,...,bTO] and let B' = [d, bg, . . . , b^]- Without loss of 
generality, we assume B' has rank to, otherwise we can simply reorder the 
basis vectors. Let A e Z"^" such that B' = BA. We now decompose A into 
Hermite normal form, which can be done in polynomial time in the input size 
and the dimension [21]. That is, wc find a unimodular matrix U € Z"^" and 
an upper triangular matrix T E j/i-Km g^,}^ ihdX A = UT; therefore, B' = 
{BU)T. There are several algorithms to compute Hermite normal form. For 
a worst case complexity bound, wc use Storjohann and Labahn [34] with run 
time 0{n'^M{nL)) where L is a bound on the maximum binary encoding length 
of each entry of A, and M{t) is the time required multiply two numbers of 
size [t] . The entries of A arc all I's, O's, and A^'s. Since U is unimodular, 
BU = [bi, . . . , hfn] is a basis for A. Since T is upper triangular, we find that 
Tiibi = d, and because d is primitive and bi G A, we have bi = d. Thus 
{d, b2, • • • , bm} is a basis for A. □ 

For the case where A = Z" , we can choose B = [ei , . . . , e„] . For Lenstra's 
algorithm, ||d||2 < n, hence the time complexity simply becomes n'-^^^\ 
After choosing a flatness direction d, if the width of the inner ellipsoid is 
smaller than oj{n), wc will project into hyperplanes perpendicular to d. Ac- 
cording to Lemma 3.8, we compute vectors bi, . . . , b„_i G Z" such that B = 
[bi, . . . ,b„_i,d] is a basis for the lattice Z". We now consider the projection 
of Y into the hyperplane d^x = t, 

Ft = {xeM"-i ■.B[x,t]eY}. 

This step must ensure the set Yj is also a convex set in class C to allow for 
the recursion to work. Since d is a flatness direction of E{A, a), we know that 
]d^(x — a)| < uj{n)(5{n) for every x e E{A, a) and it suflaces to choose t in the 
set 

{z + [d'^aj : \z\ < oj{n)P + l,zG Z}, 

which has fewer than 2u{n)^ + 3 elements. This means that if the algorithm 
runs to its full extent, the total number of subcases it will have to evaluate is 

n 

worst # of subcases = J|(2w(i)^(i) + 3). (3) 

Heinz and Khachi^an and Porkolab follow Lenstra's original algorithm which 
uses u){n) = 2^^^ ^ and /3(n) = 0{n^^'^), which leads to a total number of 
subcases 

n 

original worst # of subcases = Y[{2"^'">0{i-^^''-) + 3) = 2''("'). 

i=l 
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Our new algorithm has three important features. First, it applies the SVP 
algorithm of Micciancio and Voulgaris to obtain a flatness direction in n'-^*^^^2^" 
time as opposed to the previous 2'^("'°s"^ time of Kannan. We use the best 
known flatness constant for ellipsoids, Lu{n) = n. And we achieve an 0{n)- 
rounding by using m := n2" test points for our separation oracle. This is the 
first time, to our knowledge, that this choice has been made, and we point out 
the important fact that making this choice improves the exponential coefficient 
in the final complexity. In our algorithm, the worst case number of subcases 
reduces to 

n 

modern # of subcases = [|(2(i)(0(i)) + 3) = 2"{n\f < 22"i°g2(")+". (4) 

i=l 

4 Lenstra-type Algorithm 

Here we state a modern Lenstra-type algorithm for a class C of convex sets and 
we comment on the overall complexity. This algorithm requires that C be closed 
under projection into afHne hyperplanes in M". 



Input: A convex set Y C K" in class C. 

Output: A point x* S F n Z" or confirmation that no such point exists. 
PROCEDURE: 

1. Bounds: Determine R G Z_)_ and an e > such that Y C B{0,R) and if 
F n Z" ^ 0, then vol(y) > e. 

2. Ellipsoid Rounding: Compute an ellipsoid £'(A, a) for such an e such 

that either 

(a) Y C E{A,r) and vol(£(A,a)) < e, or 

(b) E{A, a) is a /3-rounding of Y. 

If we are in case (a), then no such point exists. 
Otherwise we proceed as we are in case (b). 

3. Flatness Direction: Compute a flatness direction d e Z" of E{-^A,a.). 
Either 

(a) w{E(-w2A,a)) > Lo{n), then there exists a point x* G EijjA.a) C 
F n Z", which we can compute by solving a closest vector problem, 
or 

(b) proceed knowing that 'w{E{A,a.)) < Lo{n)P{n). 

4. Sublattice: Compute vectors bi, . . . , b„_i e Z" such that 
B = [bi, . . . , b„_i, d] is a lattice basis for Z". 

5. Project: For each t € {2 + [d^^aj : \z\ < uj{n)/3{n) + 1, z G Z}, solve the 
n — 1 dimensional integer feasibility subproblem on the set 

Yt = {xe : B[x,t] e Y}. 
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Considering Lenstra's algorithm in the form presented here, we find it has 
time complexity of 

(n + {R) + (e))'^(^>^ (Shallow Cut + m x feasible test) n^^^j 2^"^ 2^nios^{n)+n ^ 

Ellipsoid Rounding Shallow Cut Oracle SVP/CVP Sublattice Subcases 

where we have left the shallow cut and the feasibility test as unknowns since 
they are dependent on the class C. 

5 Quasiconvex Shallow Cuts 

In this section we will show that the modern Lenstra-type algorithm can be 
applied to convex sets given by quasiconvex polynomial inequalities. We will 
begin with a contribution on efficiently encoding polynomials to exploit sparsity. 
We will then review properties of quasiconvex polynomials that will be useful 
for making shallow cuts and present our shallow cut oracle. 

5.1 Polynomial Encoding 

In this paper, we allow our complexity results to vary based on the encoding 
scheme chosen for the polynomials. Multi- variable polynomials can be presented 
in a list of the coefficients of all the monomials up to degree d, requiring a large 
storage space. This is typically referred to as a dense encoding. Under this 
scheme, the following remark holds. 

Remark 5.1 (Remark 2.1 in [15]). Let F e Z[x] be a polynomial of total degree 
d at most with integer coefficients of binary length bounded by /. Moreover, 
let X G Q" be a fixed point with binary encoding size (x) < r. Then there is 
an algorithm with time complexity (Irn)'^^^^ d'-"^^^^ and output-complexity (Z + 
r)((in)'^(^) which computes the value of the function F and the gradient V-F at 
the point x G Q". 

This time-complexity, however, is too pessimistic; for example, it seems to re- 
quire n*^("^ time to evaluate a monomial of degree d = n. 

An alternative is sparse encoding, where monomials are listed with their non- 
zero exponents and their coefficients, allowing for a more concise representation 
for short polynomials and a more refined time-complexity analysis. Polynomi- 
als and their gradients can then be evaluated in (IrdMn)'-'^^^ time, where M 
is a bound on the number of monomials in the polynomial. This scheme is 
potentially problematic in Lenstra's algorithm because each subproblem is re- 
alized by intersecting our region with a hyperplane, which would cause a loss 
of sparsity (fill-in). For instance, if the given polynomial is and our hyper- 
plane is a;„ = a;i -|- • • • -|- Xn-i + 1, then in the reduced dimension it becomes 



16 




(a) -1 1 (b) -1 1 



Figure 8: (a) Convex polynomial f{x) — x^, (b) quasiconvex polynomial g{x) = 

(xi + • • • + Xn-1 + I)''- We note that in the algorithm, we never expand these 
expressions, allowing sparse encoding to continue to be useful. We instead leave 
the polynomials alone and store coordinate transformation matrices at each step 
and then compute the coordinates in the original space to input into the polyno- 
mials. Gradients are computed via the chain rule. For context, this is discussed 
in more detail in Remark 5.6. 



5.2 Quasiconvex Polynomials 

A function F: M" — > M is called quasiconvex if all the lower level sets {x G 
M" : F(x) < a},a e M are convex subsets of M". Although quasiconvex func- 
tions are not necessarily convex, all convex functions are quasiconvex. We fol- 
low [15] for a review on quasiconvex polynomials. 

Lemma 5.2 (Section 4.1, Remark 1 in [7]). Let F e ]R[x] he a quasiconvex 
polynomial, a e M" a fixed point and b G M", b 7^ a fixed vector. If the 
polynomial i^(a-l-Ab) in A G M is strongly decreasing (or constant, respectively), 
then F{ii.-\- Ab) is strongly decreasing (or constant, respectively) for all x G M". 

This lemma does not necessarily hold if the function is not a polynomial. 
Consider an example from [9], /(x) = max{i: Xi 7^ 0}. This is quasiconvex 
because all the lower level sets are subspaces, for example {x: /(x) < n} = 
{x: Xn = 0}. This is a counterexample since /((xi, 1)) = n (a con- 
stant), but /((xi, x„_i, 0)) can vary with the remaining input. 

Lemma 5.2 can be used to determine if a quasiconvex polynomial is constant. 

Corollary 5.3 (Corollary 2.3 in [15]). Let F G M[x] be a quasiconvex polynomial 
of degree d at most, a G M" a point, and let the set {bi, . . . ,b„} C M" he a 
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basis ofW"'. If for every i = 1, . . . ,n, there are pairwise distinct real numbers 
Ail, . . . , Xid e M satisfying VF(a + Xijhi) = for all j = 1, . . . , d, then the 
polynomial F is constant. 

The following lemma is important for generating shallow cuts. 

Lemma 5.4 (Lemma 2.4 in [15]). Let F e M[x] be a quasiconvex polynomial 
and let y e M" be a fixed point. If F{y) > and VF{y) ^ 0, for every other 
X e M" that satisfies F(x) < 0, we have that 

VF{y) ■ X < VF(y) • y. 

As mentioned before, the class of convex sets used must be closed under 

intersections with affinc hypcrplancs. We will also require that an ellipsoid 
bound Aqx < R reduce to similar ellipsoid bound. 

Remark 5.5 (Within the proof of Theorem 4.2 in [15]). Let Fq, . . . , e Z[x] 
be quasiconvex polynomials, iieZ, j4oeZ"^"a positive definite matrix, and 
Fj+i G Z[x] a polynomial defined by F,+i(x) := —R + x'^AqX, for x G M". 
Moreover, let the binary length of the coefficients be bounded by I, let d be an 
upper bound for the degree of the polynomials. Let B e Z"^" be nonsingular, 
f e Z, with entries of B and t of binary length at most l{dn)^^^\ Let 

y = {x e M" : Fi(x) < 0, i = 1, . . . , s + 1} 

and let 

Yt := {iG M""^ : B[i,t] G Y}. 

Consider the set Yt and the new coordinates xi, . . . , x„_i induced by x = B[x, t], 
fixing the last coordinate in — t and rewriting the quasiconvex polynomials in 
terms of the new coordinates. The maximum binary length of all coefficients 
belonging to the new polynomials Fq, . . . , Fg+i G Z[ii, . . . ,x„_i], is l{dn)^^^\ 
Furthermore, all new polynomials are quasiconvex since the transformation is 
linear and F^+i preserves its form for a new suitable Aq. The degree bound d and 
the number s of polynomials remain unchanged, but the number of coordinates 
reduces by one. 

Remark 5.6. Following the notation of Remark 5.5, wc will explain here how we 
evaluate the polynomials and their gradients under the sparse encoding scheme. 
Suppose k+1 of such coordinate transformations B", . . . , e Z"^" are done 
to produce the variable x"~'^ G Z""*^"-*^. Each is a block diagonal matrix 

where the last block is an identity matrix of size i. In each transformation 
we are restricting the last variable to be tj. A polynomial F transformed into 
the new coordinates we will denote as F"~'°. For a given x"~'^ G Z"~*, we can 
compute pi^-k^j^n-k-^ 

pn-k^^n-k-^ = F(x) 

where 



x = B"B"-i---S"- 
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With the product (7""*^ = i?" • • • ^"-'^ computed ahead of time, a depth first 
search allows us to store at most n of these products at any given time. The 
partial derivatives of then have a simple representation as 

f)pn—k a„ 

where C""'' is the i**^ column of C""*^. 

5.3 Shallow Cuts 

The main result of [15] is derived from Heinz's shallow cut separation oracle 
for quasiconvex polynomials. The following is an adaptation of Heinz's proof to 
allow for stronger ellipsoid roundings. Specifically, we show that his calculation 
for a basis direction to admit a shallow cut generalizes to having any direction 
admit a shallow cut. Recall that we are solving the feasibility problem over the 
set 

y = {x e K" : Ji(x) < for i = 0, 1, . . . , s} 

where all the F[s are quasiconvex polynomials with integer coefficients. Consider 
an ellipsoid E{A,a) and let {bi, . . . ,b„} C M" be an orthogonal basis of M" 
according to the matrix A (where the inner product is given by (x, y) a = x-^Ay). 
Define the afRne map r : K" ^ M" such that 

T(x):=B^(x-a) (5) 



^^=(m^'---'^^V)eK""". (6) 



where 

^■=(- 

|bi|U l|b„||A 
Thus x^x < 1 if and only if r-i(x) e E{A,ai). 

Theorem 5.7 (Shallow Cut Separation Oracle). Let c > 1 and Zei m: N — > N, 
such that 2n < m{n) < c". Then there exists a function /3: N — >■ M where 

3/2 

/3(n) = Q( iog"TO/„) ) (ind an algorithm with the following input: 

(11) sparsely encoded quasiconvex polynomials Fo,...,Fs+i € Z[x] of total 
degree d, at most M monomials in each, and whose coefficients ' binary encoding 
lengths are bounded by I, 

(12) an ellipsoid _E(y4., a) containing Y as defined in (2), where the binary 
encoding length of the columns of A and of a are bounded by r, 

and outputs one of the following answers: 

1. confirmation that the ellipsoid E{A,a.) is a j3-rounding ofY , or 

2. a vector c € Q", c ^ 0, with the property 



FC |xe]R":c^x<c'^a+^-||c|u!>. (7) 

n + 1 J 
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This algorithm runs in time- complexity s{lnrmM)^^^^ and with 
output- complexity {I + r){dn)'-'^^K 

Proof. First compute an orthogonal basis {bi,...,b„} according to A. Let 
cr > 1 (Heinz used a = 3/2). Next construct a polytope approximating 
according to Theorem 2.3 using m vertices and let y c Z" denote the set of 
non-normalized vertices. For every v € V, define 



bv:= 



|v||2||bi| 



(n + a) 

Since we cannot compute bv and Xv exactly due to square roots from the 
norms, we will approximate the square roots. This can be done with any 
root finding technique rising fixed point arithmetic, for example, the Ncwton- 
Raphson method has quadratic convergence and will find the desired approxi- 
mation within polynomial time. For a reference on numerical methods, see [10]. 
We note that 

||v||2||bi|U = V^^^JhfAhi = ^Jv^vhfAhi. 

Thus we will approximate the reciprocal of that square root within an accuracy 
of 5 = e/(||^^/^||oo||bmax||2\/^) where ||^||c» is the maximum row sum of A 
and bmax = ai'gmax{||bi||2 : i = Let < (Sv,i < S he the exact 

error on each approximation. Let bv & Q" and Xv € Q" denote the rational 
approximations of bv and Xv respectively. 

Since V is symmetric across all axes, it suffices to store only the vertices in 
the first orthant. This exponentially reduces the number of root approximations 
necessary. 

Case 1: (fi-rounding) Suppose Xv € 5^ for all w & V. We will show that E{A, a) 
is a /3-rounding of Y. 

First observe that 

T(a + ylbv) = B^Ah^ 

hjAh, 



^ ht Ah, ( 1 , ^ 

~hC'\M\A Vl|v||2||b,|U+'-' 
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where e K" is the i^^ unit vector. Hence, we have 

V 



V 2 



T(a + Abv) 



i\\A 



= max Vi^v-JIb 

i— l,...,n 

< max 5v,i||bj||^ 

z— l,...,n ' 
<<5||Al/2||^||b„a.||2 

< e/^/n. 



Since II • lloo < f/v^ implies || • II2 <i we see that {S'^^bv : v G is an 
e-approximation of V (after normalizing to the unit sphere). Therefore by Corol- 
lary 2.5, there exists a. f3 = 0(-\/n/log(m/n)) such that conv(ii') C 0). 

Letting /3 = /3(n + ct), it follows that 

E{j,A,a) = T-\E{j^I,0)) C conv({iv : v e V}) C Y C E{A,ai). 

Case 2: (Shallow cut) Suppose that Xv ^ ^ for some v e V. We will show 
that there exists the desired hyperplane. 

For some F € {Fq, . . . , F^}, we know that i^(xv) > 0. 

Case 2.1: F{a.) < 0. We will find a point near Xv with non-zero gradient. 
Pick scalars Ai , . . . , A^, such that 

n+1 , ,1 

— — < Ai < ■■■ < Ad < . 8) 

n + a llbvlU 

Since the inequalities F{a) < and F{5c^) > are valid, the polynomial 

F{a. + ^^^Abv) is of degree d at most and not constant with respect to A. 
Therefore, for some 1 < fc < d, we may choose a point y e M" satisfying 

y = a + ^— AfeAbv and VF(y) ^ 0. 
n+1 

Define c :— VF(y)"^. Note that y ^ Y since Xv is a convex combination of a 
and y. Thus for all x e F, 

c^x < c^y = c^a+^— AfcC^^bv < c^a+^— ^.=^^ < c^a+^— ||c||a. 

n+1 n+1 llbvlU n + 1 

The last inequality comes from the Cauchy-Schwarz inequality for the scalar 
product generated by the matrix A. 

Case 2.2: F{a) > (i.e., a ^ Y). We can now completely ignore the fact 
that Xv ^ F and instead just find a unit direction point with non-zero gradient. 
If VF(a) 7^ 0, then we can simply use c = V-F(a). Otherwise, we may need to 
examine up to nd points not in Y. 

For this reason, pick scalars A^i, . . . , Aj^ for each i = 1, . . . ,n such that 

< Ail < • • • < Aid < jwrrr 
bi U 
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once again using the Ncwton-Raphson method to approximate the roots, but 
this time the same precision is not required. 
Define the set 

C = {a ± jXijAhi : j = l,...,n} , 

where < 7 < -j^. By Lemma 5.3, if V^'(y) = for all y € C, then F is 
constant. If so, any point c e Q" will suffice to output. 

Otherwise, we will show that there exists a point in C \ F with non-zero 
gradient. 

We remark here that we only need to choose 7 = to generate a shallow 

cut, but choosing a smaller 7 will allow for a deeper cut and therefore a faster 
convergence of the shallow cut ellipsoid method. Unfortunately, this does not 
improve the theoretical complexity of the algorithm. See [13, chapter 4] for 

more about the complexity of various cuts. 
For every i = 1, . . . , n, define the finite subsets 

C+ := {a + 7Ay Abj : j = 1, . . . , n} , 
C~ := {a - jXijAhi : j = 1, . . . ,n} . 

For every i = 1, . . . ,n, since a is a convex combination of points in and C~ , 
at least one of the sets (lY 01 C~ fl F is empty. Thus there exists a point 
y G C such that, similar to Case 2.1, 

F{y) > and Vi^(y) ^ 0. 

Define c = VF(y). The remaining calculation is similar to the one above. □ 

This proof differs from Heinz's in three ways. Most obviously, we generalize 
the test points for a rounding based on any chosen m > 2n. Secondly, we 
generalize Heinz's cuts to show that any point y with r(y) € i?(^^,0) and 
'^P{y) 7^ generates a shallow cut. Lastly, in case 2.2, Heinz's proof requires 
Xij bounded from below as the previous Aj coefficients were. As we show, that is 
not necessary, and is even disadvantageous. Allowing a different 7 for a deeper 
cut is a new idea. 

Corollary 5.8. Let c > 1 and let m: N ^ N, such that 2n < m.{n) < c". Then 
there exists a function /3: N — >■ M where (3{n) = Q( iog"^/„) ) "^(^ algorithm 
with the following input: 

(11) sparsely encoded quasiconvex polynomials Fq, . . . , Fg e Z[x] of total de- 
gree d and of at most M monomials, 

(12) an integer R and a positive definite matrix Aq e ^nxn define 
Fs+i(x) —R + x^Aqx. Let I be a bound on the maximum binary encod- 
ing length of the coefficients of Fq, . . . , i^s+i; 

(13) a positive number e e Q, 

which outputs a positive definite matrix A e Q"^" and a point a e Q" such that 
one of the following holds: 
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Figure 9: Case 2.1: If -F(a) < and -F(xv) > 0, then F restricted to the Une 
through those two points is not constant. One of the test points on that hne 
must have a non-zero gradient and not he within Y . 




Figure 10: Case 2.2: If F{a.) > we search for a point near a in a unit direction 
that has gradient non-zero and is not in Y. 



23 



1. Y C E{A, a) and vo\{E{A, a)) < e, or 

2. the ellipsoid E{A,a) is a ^-rounding ofY. 

This algorithm runs in time- complexity s{lnmdM {e))^^^^ and with 
output- complexity {I + {e)){dn)'^^^\ 

Proof. The proof is similar to [15, Corollary 3.4]. □ 

6 Lenstra-type Algorithm 

Now that we have a shallow-cut separation oracle, our Lenstra-type algorithm is 
almost ready to be applied to quasiconvex integer minimization. Before arriving 
at the proof of main theorem, we state the following lemmas related to the 
convergence of the shallow-cut ellipsoid method and conclude that we can solve 
the feasibility problem. 

Lemma 6.1 (Lemma 4.1 in [15]). Let Fq,...,Fs+i G Z[x] be polynomials, 
R G Z an integer, Aq G Z"^" a positive definite matrix, and Fg+i € Z[x] a 
polynomial defined by _F5+i(x) := —R + x^^AgX for x G K". Moreover, let 
the binary encoding length of the coefficients be bounded above by I, let d be 
an upper bound for the degree of the polynomials, and let the set Y contain an 
integer point x G Z". Then there is a positive rational number e G Q which 
bounds the volume < e < vol{Y) such that its binary length (e) is in the class 
l{dnf(^l 

Lemma 6.2 (p. 27 in [6]). // a minimum point exists to problem (1), then there 
exists a ball of radius R^ G Z containing such a point, where the binary length 
ofR^ is Id'^^'^l 

The following theorem solves the feasibility problem. 

Theorem 6.3. Let Fo,...,Fs G Z[x] be quasiconvex polynomials, R > an 
integer, Aq G Z"^" a positive definite matrix, and -Fs+i G Z[x] a polynomial 
defined by Fs+i{x.) = —R + x^Aqx for x G R". Let d be an upper bound for 
the degree of the polynomials Fq, . . . , Fs+i, presented as a sparse list of m,ono- 
mials with at most M monomials, and let the binary length of the coefficients 
be bounded by I. Moreover, consider the set 

y :={xe]R" :Fi(x) < 0, i = 0, 1, . . . , s + 1} (9) 

Then there is an algorithm with time- complexity s(dM^)'-'(^)2^"'°S2(")+'^("') which 
computes a point x* e F fl Z" or confirms that no such point exists. 

Proof. This is a direct result of the complexity analysis done in section 4 in 
conjunction with Lemma 6.1, Lemma 6.2, and Corollary 5.8. □ 

From the approach of this method, due to the recursion, it is unlikely possible 
to obtain a better time-complexity than 2*^^" in terms of the dimension. 
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An open problem is then to find an algorithm that is single exponential in the 
dimension. 

We now provide an outline for the proof of Theorem 1.1, which follows from the 
same reasoning as [15, Theorem 5.1]. 

Proof of Theorem 1.1. Again, from [6, p. 27] that if a minimum point x* exists, 
then there exists a ball of radius i?* S Z containing such a point, where the 
binary length of i?* is r = ldP^"\ Therefore, F{x* ) has binary length bounded 
by («rdMn)°(i) = (/Mn)°(i)d'^("). We define the polynomial F,+i e Z[x] by 
Fs+i{x): = — i?2 + x'^x for x e M", making the problem bounded. To solve 
the optimization problem, we compute the smallest integer z* € 1^ such that 

{x e Z" : Fo(x) - z* < 0, and F,(x) < for i = 1, . . . , s + 1} 7^ 0. 

We then apply binary search on z * within binary length of (?Mn)°(i)d°("), 
testing for integer points using Theorem 6.3, to find an optimal z* and obtain 
the desired time-complexity because of the bound given by [6] . □ 
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